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ABSTRACT 

We present new spectroscopic observations of the stellar cluster population of region B in the prototype star- 
burst galaxy M82 obtained with the Gillett Gemini-North 8.1 -metre telescope. By coupling the spectroscopy 
with f/Sy/ photometry acquired with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope 
(HST), we derive ages, extinctions and radial velocities for seven young massive clusters (YMCs) in region B. 
We find the clusters to have ages between 70 and 200 Myr and velocities in the range 230 to 350 km s^^, while 
extinctions Ay vary between 1-2.5 mag. We also find evidence of differential extinction across the faces of 
some clusters which hinders the photometric determination of ages and extinctions in these cases. The cluster 
radial velocities indicate that the clusters are located at different depths within the disk, and are on regular 
disk orbits. Our results overall contradict the findings of previous studies, where region B was thought to be 
a bound region populated by intermediate-age clusters that formed in an independent, offset starburst episode 
that commenced 600Myr-l Gyr ago. Our findings instead suggest that region B is optically bright because of 
low extinction patches, and this allows us to view the cluster population of the inner M82 disk, which probably 
formed as a result of the last encounter with M8 1 . This study forms part of a series of papers aimed at studying 
the cluster population of M82 using deep optical spectroscopy and multi-band photometry. 
Subject headings: galaxies: evolution — galaxies: individual (M82) — galaxies: photometry — galaxies: spec- 
troscopy — galaxies: starburst — galaxies: star clusters 



1. INTRODUCTION 

The extensively studied galaxy M82 is a local example of 
a nuclear starburst galaxy. O'Connell & Mangano (1978, 
hereafter OM78) first catalogued the complex star-forming 
regions seen in ground-based images of the M82 disk, and 
introduced the nomenclature A-H. Region B is the brightest 
region in the disk and is located 350-1050pc north-east of the 
nucleus. OM78 and Marcum & O'Connell (1996) find that 
the integrated spectrum of this region is indicative of a fos- 
sil starburst region as it shows the characteristic 'E+A' post- 
starburst spectrum, suggestive of a truncated burst of star for- 
mation occurring 100- 1000 Myr ago. Moreover, they find that 
the intrinsic brightness of region B is such that the burst must 
have been of comparable intensity to the present starburst in 
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the nucleus. 

More recently, de Grijs, O'Connell, & Gallagher (2001) 
studied the M82-B cluster population using photometry ob- 
tained with the Wide Field and Planetary Camera 2 (WFPC2) 
and the Near-Infrared Camera and Multi-Object Spectrome- 
ter (NICMOS) onboard the Hubble Space Telescope (HST). 
They identified 113 clusters and, by estimating ages and ex- 
tinctions from BVI photometry, they found that a concen- 
trated episode of cluster formation occurred 400- 1000 Myr 
ago with a peak at 600 Myr Subsequent to this analysis, 
de Grijs, Bastian, & Lamers (2003a) re-derived ages and 
extinctions by combining BVI and JH photometry; they 
find that the peak of cluster formation occurred at 1.10 Gyr 
with an age spread of 500-1500Myr. de Grijs, Bastian, 
& Lamers (2003b) used these new ages to derive the clus- 
ter luminosity function (CLF) for a fiducial age of 1.0 Gyr 
and find that it has a Gaussian shape, similar to globular clus- 
ters (GCs), rather than the standard power-law CLF found for 
populations of young massive clusters (YMCs) (e. g. Larsen 
2004; Gieles et al. 2006a). This result is surprising but sup- 
ports theoretical models advocating that an initial power-law 
distribution of cluster masses will be transformed into a Gaus- 
sian distribution because low mass clusters will be preferen- 
tially disrupted (e.g. Fall & Zhang 2001; Vesperini & Zepf 
2003; Gieles et al. 2006b). de Grijs et al. (2003b) therefore 
suggest that because M82-B is of intermediate age, it provides 
the 'missing evolutionary link' between the power-law CLFs 
found for YMCs and the Gaussian CLFs of globular clusters. 
However, for an environment as turbulent as that of M82, one 
would expect a timescale of several Gyr for preferential dis- 
ruption to produce this effect (e. g. Goudfrooij et al. 2004, 
2007, where such a timescale is calculated forNGC 1316 and 
NGC 3610 respectively). 

To date, spectroscopy of the region B cluster population 
has been extremely limited; the only published spectra are 
those of Smith et al. (2006) for the two brightest members. 
The spectra were obtained with the Space Telescope Imag- 
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ing Spectrograph (STIS) and, although they are of low qual- 
ity, the derived ages of 350±100Myr are much lower than 
the photometrically-based ages of ~0.7-6 Gyr (Smith et al. 
2006; de Grijs et al. 2003a) for the same clusters. This dis- 
crepancy hints at the possibility that region B may be younger 
and that the ages derived from the photometry are too high. It 
is important to determine accurate ages for the M82-B cluster 
population to verify its unusual CLF, and for understanding 
the cluster formation history of M82 and its relationship to 
encounters with its close neighbour M81. 

We therefore acquired new spectroscopic and photometric 
data for the M82-B cluster population. By using both tech- 
niques we are able to obtain information for a large number 
of clusters and also overcome the degeneracy between age and 
extinction that presents a hurdle in the analysis of photometric 
data (Trancho et al. 2007a,b). In this paper, we present opti- 
cal spectroscopy for seven clusters obtained with the Gillett 
Gemini 8.1-m telescope, and in a companion paper (Smith et 
al. 2007), we present new HST imaging, including ?7-band 
photometry of the clusters. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spectroscopy 

The spectroscopic data presented in this paper were ob- 
tained with the Gillett Gemini North 8.1 -metre telescope us- 
ing the Gemini Multi-Object Spectrograph (GMOS-N), as 
part of program GN-2006A-Q-38 (PI Smith). The target se- 
lection was made based on photometric criteria with the po- 
tential YMCs identified on GMOS pre-imaging, taken as part 
of the observing program. A multi-sUt mask was designed, 
using a slit width of 0'.'75, with a variable slit length that we 
tried to keep to a minimum of 6'.'5, in order to include an area 
of sky that would ensure a good background subtraction. In 
some cases however, a smaller sUt length was used so as to 
include a greater number of sources on the mask, resulting in 
a total of 39 objects being targeted across the M82 galactic 
disk. In this paper we focus on seven of these 39 clusters that 
are situated within region B. Fig. 1 indicates the positions of 
the clusters and the sUts, and the fully reduced rectified spec- 
tra are presented in Fig. 2. A list of cluster identifiers and 
coordinates can be found in Table 1 . 

The data were acquired on the night of 2006 April 5, under 
good seeing'' conditions (0'.'8 at 5000 A). We used the B600 
grating and the data were read-out in 2x2 biiming mode, re- 
sulting in a dispersion of 0.9 Apx~^. As this is multi-object 
spectroscopy, the central wavelength varies slightly with ev- 
ery slit, so the resulting spectral range of the data is ap- 
proximately 3700-6500 A. While we used all three GMOS-N 
CCDs to achieve this broad range, the sensitivity of the silver- 
mirrored chips is lower at the blue end, resulting in a poor 
S/N and loss of accuracy in the wavelength calibration below 
4000 A, thus setting a lower hmit to our effective wavelength 
range. 

The data were acquired in eight 1800 second exposures, 
totalling four hours of integration time. More specifically, 
there were two sets of four exposures that differed in terms 
of central wavelength. This is standard GMOS procedure 
which allows the observer to compensate for the inter-chip 
gaps through co-adding/merging the spectra; the two Accntrai 
were 5080 and 5 120 A. CuAr (Copper- Argon) arc exposures 

^ This is not the seeing, but the 'image quality' parameter adopted by Gem- 
ini, which incorporates both wind and telescope performance effects. 



and Quartz-Halogen flat fields were taken in between target 
exposures. Bias frames were taken as part of the Gemini base- 
line caUbrations (GCAL). 

The data were subsequently reduced using mainly standard 
IRAF^ reduction tools, as well as tasks from the purpose- 
developed Gemini-IRAF package (vl.9.0). The data were 
bias-subtracted, flat-fielded and then combined into a single 
image, while correcting for discrepancies in quantum effi- 
ciency between the three detectors using the QE Gemini IRAF 
tool. The wavelength calibration was performed using the ob- 
tained CuAr arc frames, resulting in residuals typically less 
than 0.15 A. 

One comphcation that arises when observing in multi-sht 
mode without an active atmospheric diffraction correction de- 
vice is differential refraction; this causes wavelength depen- 
dent slit-losses, as it is not possible to keep the slits aligned 
with the parallactic angle for the duration of an exposure. To 
alleviate this effect, the combined data were corrected for par- 
allactic angle using a purpose-built IDL (Interactive Data Lan- 
guage) routine written by B. W. Miller (2006, priv. comm.), 
based on Filippenko (1982), where the effects of atmospheric 
differential refraction are discussed. 

At this point we extracted all 39 spectra from each multi- 
dimensional MOS frame and merged all available exposures 
of each target. In most cases we decided not to use all avail- 
able spectra, as some exposures were of poor quality (caused 
by seeing degradation) and thus lowered the quality of the 
merged spectrum. 

Finally, the spectra were flux-calibrated using a response 
function derived from the standard star Wolf 1346. It should 
be noted here that the observations used were taken on a dif- 
ferent night (2006 May 25), due to problems that arose with 
the exposure of the original flux standard, Feige 34. This un- 
fortunately led to an increased uncertainty in the absolute flux 
scale. 

Having completed the extraction procedure, we measured 
the width of the sky-line at A5577, in order to measure 
the spectral resolution of our data, which we found to be 
3.5 A (FWHM). We also tested the wavelength calibration by 
measuring the exact position of this line and found it to be 
well within the predicted errors. The signal-to-noise ratios 
(S/N) of the spectra were determined using two continuum 
windows between 4400-4800A and 5000-5800 A. We found 
values ranging from 16^4; the individual values for each 
cluster are presented in Table 3. 

2.2. Photometry 

We used imaging from the recent Hubble Heritage ACS 
Mosaic of M82 (for full details of the observations see Mutch- 
ler et al. 2007). This is a composite of six individual pointings, 
resulting in a combined field of view of 12288 x 12288 pix- 
els, or approximately 9x9 arcininutes squared, given the 0.05 
arcsec per pixel scale of the ACS-WFC. The images cover 
the entire galaxy in four bands, of which we use the F435W, 
F555W, and F814W that correspond roughly to the Johnson 
B, V and /. We performed aperture photometry of the sources 
in region B using an aperture of 10 pixels (0'.'5) and a back- 
ground annulus with irmer radius and width of 1 1 and 1 pixels 
respectively. 

We also used f/-band (F330W) ACS-HRC imaging of re- 

' IRAF is distributed by the National Optical Astronomical Observatories, 
which are operated by the Association of Universities for Research in Astron- 
omy, Inc. under contract with the National Science Foundation 
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Fig. 1 . — Finding chart for all the clusters in region B discussed in this paper. The inset shows the location of the region within the galaxy. East is towards 
the left and North is to the top of both images. The rectangles represent the slits used to acquire the data, to scale, with identifiers placed above the clusters. 
We choose to study these seven clusters because they are bright and isolated. Note the intricate structure of the gas and dust in this region. Both images are 
composites of four filters, B, V, I and Ha (taken from the Hubble Heritage HST ACS mosaic; see Mutchler et al. 2007). 



TABLE 1 

Positional information for the cluster sample. 



#" 


Alt*'' 


Previous 


RA'' DEe 


91 


34 


B2-12, H', B1-2E 


9 55 54.60 69 41 01.7 


97 


25 


B2-26 


9 55 55.65 69 41 05.9 


103 


20 


B2-41 


9 55 57.04 69 41 07.4 


108 




B2-49 


9 55 58.41 69 41 08.7 


113 


14 


Bl-6 


9 56 00.63 69 41 02.5 


126 


6 


Bl-20 


9 56 02.52 69 41 08.1 


131 


1 


Bl-28, B2-1'' 


9 56 03.40 69 41 12.3 



"Cluster identifier as in our Gemini program 
''Identifier as in Smith et al. (2007). 

'^Identification by de Grijs et al. (2001, unless otherwise stated), with 'Bl' 
and 'B2' referring to the part of the region to which the cluster belongs. 
''•'^12000 coordinates as measured on the F555W ACS mosaic image; the units 
are (h m s) for a and (° / //) for 5. 

fNotation follows OM78. 
s.'' After STIS study of Smith et al. (2006). 



gion B, obtained as part of HST program #10853 (PI Smith), 
which extends our photometric coverage to the near-UV. The 
combination of optical and UV data enables us to set a strong 
constraint on age and extinction determinations of young stel- 
lar clusters (e.g. Anders et al. 2004). The HRC data were pro- 
cessed with the standard automatic HST pipeline and drizzled 
(using the MultiDrizzle package - Koekemoer et al. 2002) to 
a scale of 0'.'025 per pixel. For a more detailed discussion 
of the data we refer the reader to Smith et al. (2007). We also 



performed aperture photometry on the HRC images with aper- 
tures of the same physical size (Of.'S as projected onto M82), 
that is, an aperture radius of 20 pixels and a background annu- 
lus of inner radius and width of 22 and 2 pixels respectively. 

For both the WFC and HRC data, we transformed the mea- 
sured flux to the VEGAMAG system using the zero-points 
given in Sirianni et al. (2005). 

Due to the extended nature of the clusters in this study, we 
derived aperture corrections for each of the images based on 
bright, isolated clusters in our field of view (namely clusters 
#91 and 97). The corrections were found by comparing the 
cluster magnitudes as measured by applying apertures of two 
different radii (10 and 40 pixels for the WFC images and 20 
and 80 pixels for the HRC images). The applied corrections 
are -0.64, -0.35, -0.33 and -0.43 magnitudes in the U, 
B, V and / bands respectively. We note, however, that since 
the size and profile of young clusters vary (e.g. Larsen 2004; 
Bastian et al. 2005) these corrections may not apply precisely 
to all clusters studied, although we expect that the different 
cluster sizes will not affect the derived colours. 

Note that throughout the paper we use the Johnson U BVI 
notation for sake of simplicity, but we do not convert magni- 
tudes to this system. 

3. ANALYSIS AND RESULTS 

The combination of t/BV/ photometry and high S/N optical 
spectra allows us to determine the age and extinction of each 
of the seven clusters in our sample. Throughout this section 
we assume that the clusters have solar metallicity, in agree- 
ment with ISM studies of M82 (McLeod et al. 1993). Addi- 
tionally, we adopt a standard Galactic extinction law (Savage 
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4000 4500 5000 5500 6000 

Wavelength / A 

Fig. 2. — Continuum-normalised Gemini spectra of the seven clusters presented in this paper. The wavelength range is 4000-6100 A for all spectra. We use 
an arbitrary flux-scale for means of presentation, however, all analysis was performed on absolutely flux-calibrated data in the presented wavelength range. The 
locations of the lines we use as our main diagnostic tools are indicated and labelled at the top of the plot. 
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TABLE 2 

Summary of photometrically derived results. 





F330W 


F435W 


F555W 


F814W 


Mv 


Av 


Age 


Tage 


# 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(F555W) (Myr) 


(Myr) 


91 


20.12 


18.91 


17.75 


16.07 


-12.95 


2.4" 


C 




97 


18.93 


18.81 


18.28 


17.25 


-10.56 


1.1 


180 


[140, 220] 


103 


19.21 


19.08 


18.35 


16.94 


-11.36 


1.9 


90 


[60, 150] 


108 


d 


20.04 


19.21 


17.84 


-12.24 


2.5^ 


f 




113 


20.41 


20.10 


19.35 


18.01 


-10.30 


1.9 


170 


[72, 360] 


126 


19.77 


19.34 


18.64 


17.47 


-10.70 


1.5 


260 


[190, 1200] 


131 


18.37 


18.24 


17.74 


16.61 


-11.22 


1.2 


170 


[110, 190] 



"The minimum and maximum acceptable values for the cluster age accord- 
ing to the 3DEF method. 

"Value estimated from the extinction map in Fig. 7. This cluster is discussed 
in § 3.4. 

■^•'Not possible to obtain a photometric age because of differential extinction 

(see § 3.4). 
''Cluster is not visible in the (7-band. 

■^Estimated from BVI photometry assuming spectroscopically derived age. 



&]VIathis 1979). 

We note here that the assumptions on the metallicity and ex- 
tinction laws will not have large impacts on our conclusions. 
In particular, if we were to use the starburst extinction-law of 
Calzetti (1997) the photometric ages and extinctions would 
remain virtually unchanged. 

3.1. Ages and extinctions from photometry 

As a first method to determine the ages and extinctions of 
the clusters we compare their colours to simple stellar pop- 
ulation (SSP) models. We choose the GALEV SSP models 
(Anders & Fritze-v. Alvensleben 2003), with solar metallic- 
ity and a Salpeter IIVIF. These models use the Padova stellar 
isochrones, and the evolution of magnitude is given as a func- 
tion of time directly in the VEGAMAG system for HST filters. 
In that way there is no need to convert the measured clus- 
ter magnitudes to the more standard Cousins- Johnson system, 
which would introduce additional errors. 

In Fig. 3 we show the U—B vs. V—I colours of the clusters 
in region B for which we have UBVI photometry and spectra. 
Additionally, we demonstrate the evolution of the GALEV 
SSP models by labelling six ages, and we show extinction 
vectors of magnitude Ay =3.0 starting from these ages to 
guide the eye. From this plot we can deduce that all but one 
of the clusters (#91) have ages between 50 Myr and 1 Gyr. 
Cluster #91 (cluster H, following notation by OM78) will be 
discussed in § 3.4. 

However, as the SSP colours do not evolve monotonically 
in colour-space, the U — Bvs. V — I diagram provides sim- 
ply one of the possible colour projections of the photome- 
try. In order to fully exploit the available data, we have used 
the 3DEF method presented in Bik et al. (2003) and further 
tested in Bastian et al. (2005). Briefly, this method compares 
the photometry of each cluster to that of the SSP models, to 
which it applies an extinction (Ay) between and 5 mag, in 
steps of AAv = 0.02 mag and calculates a xt (reduced x^)- 
The model, i.e. the combination of age and extinction with 
the lowest is selected as the best fit and the errors are as- 
signed by the extrema of the set of models which satisfies the 
condition xl < X^,bcst + 1- 

The photometrically derived best fit ages range between 90 
and 260 Myr and we find Ay extinction values in the range 
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Fig. 3. — Colour-colour plot for six of the seven clusters in M82-B. The 
solid line traces the GALEV SSP model colours, with asterisks marking 
model track ages of 10, 50, 100, 250, 500 and 1000 Myr. The data-points 
represent the cluster colours and associated errors, as derived from multi- 
band photometry. The arrows represent extinction vectors, along which the 
'footprint' of the cluster may be traced, in order to estimate its age and ex- 
tinction; the length of the vectors is 3 magnitudes. The large V — I error on 
cluster #1 13 is due to contamination from a neighbouring red source. The er- 
rant cluster #91 has an underestimated t/-band flux (see Sect. 3.4) and cluster 
#108 is not included because it is not detected in the [/-band. 



TABLE 3 

Spectroscopically derived age (r) and radial velocity (u) 

measurements (HELIOCENTRIC). 





S/N" 


VR 


VR mif 


tivisrm'^ 


TCCM 




# 


(...) 


(kms-1) 




(Myr) 




91 


44 


230 ± 20 


270 ± 10 


180 


190 


[110, 250] 


97 


34 


290 ± 20 


285 ± 10 


180 


200 


[130, 320] 


103 


23 


270 ± 20 


270 ± 10 


710 


100 


[30, 180] 


108 


16 


340 ± 40 


200 ± 30 


140 


140 


[60, 250] 


113 


24 


280 ± 20 


225 ± 10 


100 


100 


[20, 180] 


126 


24 


350 ± 30 


220 ± 20 


200 


200 


[130, 280] 


131 


40 


260 ± 10 


245 ± 10 


80 


70 


[40, 160] 



"The average of calculations for two different continuum regions, A4400- 
4800 and A5000-5800. 

"The errors quoted here are adopted as a 0.25 px measurement error, which 
roughly corresponds to lOkms"^. In all but two cases the agreement be- 
tween the two individual doublet line measurements gives rise to an error 
that is lower than the adopted systematic error. 

'^•''The model-spectrum residual MSRM and cumulative CCM best fit ages 
(as explained in § 3.2). 
'The minimum and maximum acceptable ages according to the CCM. 



1.1-2.5. The values for each cluster are shown in Table 2. 

3.2. Spectroscopy 

While photometry presents a useful tool for constraining 
ages and extinctions, its limitations (such as the need for de- 
tection in four filters spanning 5000 A) show even in a sam- 
ple of seven clusters. Spectroscopy has long been recog- 
nised as providing the most accurate means of determining 
cluster ages using SSP models. We have measured ages 
for our cluster sample using the solar metallicity models of 
Gonzalez Delgado et al. (2005, hereafter GD05) in the age 
range 4 Myr-1 Gyr, which are designed specifically for young 
clusters. The parameters of the models match those assumed 
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for the photometric fitting (solar metallicity, Salpeter IMF, 
and Padova stellar isochrones), and their wavelength range of 
3000-7000A is well matched with that of our observed spec- 
tra. All model spectra have been degraded to match the ob- 
served spectral resolution of 3.5 A. 

We prefer using these theoretical models over the more es- 
tablished Bruzual & Chariot (2003) templates for a number 
of reasons. First the good age resolution means cluster ages 
can be determined with higher precision. Second, the high 
sampling resolution of 0.3 A allows us to downgrade the res- 
olution of the models, rather than that of the data, in order to 
match the two. 

We determine the ages of the clusters using two separate 
statistical methods, a model- spectrum residual method (here- 
after MSRM) and a cumulative method (CCM). We choose 
for these clusters to investigate the fit of the models on the 
lower Balmer series spectral lines. We compare the models 
and observations directly by first rectifying the spectra based 
on fits to the continuum in two windows outside of each line. 
Fig. 4 demonstrates the Balmer line fits as applied to two of 
the clusters. By fitting both H/3 and H7, we can avoid having 
to make any presumptions about the shape of the Une in cases 
where there is noise or nebular emission. We use many as- 
pects of the feature in order to determine the best fit; we check 
the agreement in line width, how well the overall line profile 
is fit and how closely the depth of the line is fit or predicted. 
In most of the clusters we find emission superposed on the ab- 
sorption line. This emission presumably comes from diffuse 
ionised gas in the region and is not centred at the same wave- 
length as the absorption, which means that it is not directly 
associated with the clusters. 

The main obstacle that arises when using automated statis- 
tical methods is the age degeneracy of Balmer line strength, 
i.e. multiple models of very different ages have similar Une 
strengths. This means that the ages cannot be determined by 
simply finding the model that best predicts the depth of the 
line. In addition to the age degeneracy, the observed spectra 
may also feature superposed emission, so the fit would not in 
such cases reflect the actual absorption line profile. 

The first method employed (MSRM) fits a model and cal- 
culates the minimised absolute value of the model-spectrum 
residual. The best fit is then chosen as the model with the 
lowest mean residual. The second method (CCM) uses all 
available data-points to calculate the cumulatively for each 
model fit. This offers a more realistic interpretation of a good 
fit to the entire line profile. In that way, the best-fitting tem- 
plate will be the one that traces the Balmer line profile, in- 
cluding the wings, thus overcoming the observed fine strength 
degeneracy. 

As mentioned above, most of the clusters contain emission 
features within the Balmer absorption lines. Therefore, in fit- 
ting the data we have carefully avoided regions that are con- 
taminated by emission. This is shown in Fig. 5, where the top 
panel shows the H7 line for cluster #97 (black solid line) and 
the two dashed boxes show where the fit was carried out. We 
would like to note here that including the profile centre in the 
fit causes degradation of the fit quality, however it leaves the 
cluster age distribution largely unaffected, with ages shifting 
by an average of ~20 Myr towards younger values. The bot- 
tom panel of Fig. 5 shows the results of the two fitting meth- 
ods applied as a function of age. The green (solid) line shows 
the CCM and the red (dashed-dotted) line shows the MSRM 
and both fits are carried out within the designated boxes. The 



MSRM is more sensitive to line strength, and has two min- 
ima for very different ages, whereas the CCM overcomes this 
obstacle in most cases and has one clear minimum. 

The error on the best fitting age is calculated in a sim- 
ilar way as for the 3DEF method, that is, based on the 
extrema of the set of models which satisfies the condi- 
tion < 2 X Xmin- ^'^^^ here that a second minimum 
does tend to appear in the lower S/N spectra (clusters #103, 
108, 113 and 126). In these cases the very bottom of the sec- 
ond minimum 'trough' is low enough to satisfy the mentioned 
error condition. Nonetheless we are able to confidently reject 
the higher age values by inspecting the resulting spectrum- 
model fit, as it shows a clear disagreement in the overall line 
profile. The derived ages and errors are given in Table 3. 

For the brightest cluster #91, it is possible to use a third 
technique to further reinforce our results, the ratio of He i, 
A4471 to Mgll, A4481; by comparing the ratio between 
spectrum and model we can get one more estimate of the 
age of this cluster. This ratio is used widely as a spectral 
type determination method in stellar astrophysics (Walbom 
1972; Walbom & Fitzpatrick 1990) and it provides an en- 
tirely extinction-independent age determination method (be- 
cause of the proximity of the two lines). In order to overcome 
sampHng-related errors that may arise when simply detect- 
ing the line minimum, we extrapolate the minimum by fit- 
ting Gaussian curves. This method has been successfully ap- 
plied to the significantly reddened cluster M82-F (Bastian et 
al. 2007). Using this technique on cluster #91, we obtain an 
age estimate of 180 Myr, which is in accord with aU other age 
measurements. 

The spectroscopic age measurements are summarised in Ta- 
ble 3. In all, we find ages in the range 80-200 Myr. We have 
adopted an age for each cluster as the average of all avail- 
able spectroscopic age determinations that provide a good fit. 
The ages are listed in Table 4. These results disagree with 
work previously published on region B (de Grijs et al. 2001, 
2003b), where it was suggested that the ages of most clusters 
in the region are in excess of 400 Myr In the de Grijs et al. 
(2003b) sample of 113 clusters, there are only four with de- 
rived ages younger than 400 Myr. This age difference and its 
implications are discussed in more detail in § 4.2. 

3.3. Radial Velocities 

The clusters observed reside in a region rich in diffuse 
ionised gas. As a result, there is a small amount of Balmer line 
emission from their environment and therefore the Balmer 
line profiles of the clusters are a blend of the deep stellar ab- 
sorption lines with narrow nebular emission lines, which in 
some cases are off-set. This makes the precise determination 
of the profile centre by inspection of Balmer lines practically 
impossible, therefore we resort to automated routines which 
also focus on other features to calculate the radial velocities 
of the clusters. To determine the velocity of the neutral gas in 
our line of sight, we use the Na I D interstellar lines. 

In order to derive the best possible estimates for the radial 
velocities we used a number of measurements (three to six 
measurements depending on the S/N), each one for a differ- 
ent spectral region and combined all results to get a value for 
the velocity (the arithmetic mean) and its error (the standard 
deviation). 

The first method employed was the FXCOR cross- 
correlating routine in IRAF, which we applied to various re- 
gions of the spectrum, each time including a different set of 
absorption lines. We also used the Penalized Pixel Fitting 
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Fig. 4. — H7 and H/3 line fitting for two of the clusters in region B. Tlie over-plotted models have been chosen to demonstrate the variation of Balmer line 
profiles with age and the effect on the fit to an observed spectrum. The top panel shows cluster #91, which stands out in Fig. 3, as it fails to meet the model track. 
We discuss the possible reasons for this in § 3.4. A first inspection of this cluster might offer an age estimate as young as 50Myr, however, care must be taken as 
there is blended, embedded, off-set emission from the surrounding gas. Comparing the resolution width to that of the instrumental sampling, it is obvious that the 
Balmer lines are well resolved. The bottom panel shows the line fits for cluster #131, which we believe to have an age of r^SOMyr. Note the slightly red-shifted, 
strong superposed emission. 
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Fig. 5. — Statistical fits for cluster #97. The top panel shows the region 
of the observed spectrum in black (wavelength in A) with the best fit mod- 
els over-plotted; the dashed green and dashed-dotted red lines correspond 
to the CCM (cumulative method) and MSRM (model- spectrum residual 
method) best fits respectively. The blue boxes represent the regions of the 
spectrum considered for the fit. The bottom panel shows the statistical mea- 
sures as a function of age; the solid green line traces the and the dashed- 
dotted red line the absolute value of the model-spectrum residual, with the 
vertical lines denoting the best fit model age for each method. In the case of 
cluster #97, this occurs at 200 and 180 Myr for the CCM and MSRM respec- 
tively (this applies to the H7 fit). The horizontal dashed green line denotes 
the Hmit of values within 2 X Xmin ^^^^ ^ ™ error estimate. The 

y-axis traces the model-spectrum residual, and the curve has been scaled 
to 1/20 of its original values in order to fit in the same plot. 



Method (pPXF), taking care not to include the emission com- 
ponent of the Balmer Une profiles in the fit. We will not pro- 
vide a detailed description of this method here; instead, we 
refer the reader to Cappellari & Emsellem (2004). In brief, 
the method simply cross-correlates selected regions of a spec- 
trum with a template in order to derive the kinematics of the 
observed object. A significant benefit of using this routine is 
that it can fit a polynomial function to the overall shape of 
the template, thus matching the effect of extinction and dif- 
ferential refraction on the observed spectra. This enables us 
to avoid making corrections or normalising the spectra before 
fitting the templates. It also means that only one measure- 
ment is required, thus it does not introduce extra degrees of 
freedom in our calculation. 

We find the two methods to provide consistent results when 
given the same templates to fit, resulting in accurate and pre- 
cise radial velocity measurements. Correcting the cluster ve- 
locities to the heliocentric frame of reference gives values 
ranging between 230 and 350kms~^ (see Table 3). This 
range corresponds to 30-150 km s^^ with respect to the cen- 
tre of the galaxy where we adopt a systemic velocity of 
200kms-i for M82 (McKeith et al. 1993). The results are 
discussed in § 4.3. 

3.4. Differential extinction on clusters #91 and #108 

In § 3.1 we briefly discussed the fact that the age of cluster 
#91 (M82-H) cannot be determined photometrically, as it lies 
beyond the area in colour-space covered by the SSP model 
evolutionary track, once extinction has been accounted for. 
This appears to be caused by an underestimation of its U- 
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Fig. 6.— ACS/HRC UBVI imaging of #91 (M82-H). The top-left, 
top-right, bottom-left and bottom-right panels represent the F330W, F435W, 
F550M and F814W band images respectively, and the orientation follows 
that of the finding chart (Fig. 1 ), with East to the left and North to the top. 
The size of each image is 1" X 1", which roughly corresponds to 18 parsecs 
at the distance of M82. The images are presented in a logarithmic intensity 
scale, while the contours show the linear scale flux levels at 20, 40, 60 and 
80% maximum intensity. Note how the U- and B-band images (F330W and 
F435W) show a highly asymmetric light profile, while in the V- and espe- 
cially the /-band (F550M and F814W) the cluster light appears to be quite 
regular. Given that the extinction will have a greater effect on shorter wave- 
length bands, we can infer that the flux in the U band is underestimated, 
hence causing the shift in colour-space that we observe in Fig. 3. 



band flux, so it also affects the analytical photometric method 
appHed (3DEF). 

In support of this hypothesis. Fig. 6 shows ACS-HRC imag- 
ing of the cluster in the U, B, V and / bands, obtained from 
the HST archive (program #10609, PI Vacca), from which 
the cluster profile appears to be highly irregular in the U- 
band, while becoming smoother and more regular in the red- 
der wavelength bands. This is the result of a complicated dif- 
ferential extinction pattern. Our measurements indicate that 
parts of the cluster are entirely obscured in the U and B bands, 
which can account for this shift in colour-space as it would 
spuriously raise the U — B value. 

In order to verify this inference, we map out the extinc- 
tion on and around the cluster pixel-by-pixel in Fig. 7; the 
colour map represents the different levels of reddening, as de- 
termined by multi-band photometry and the solid lines repre- 
sent observed isophotes in the F435W band. The technique 
we use is explained in Bastian et al. (2007). In brief, we con- 
struct a colour map by measuring the flux of each pixel in the 
B, V & I bands and therefore measure its reddening, while 
assuming the spectroscopically derived age. In principle, as 
we assume the cluster to be a SSP, we expect each pixel to 
have the same colour, therefore it is only the variation of ex- 
tinction across the face of a cluster that causes the colour to 
change. The extinction map verifies the complexity of the red- 
dening in the immediate vicinity of M82-H, showing a dust 
lane running through the eastern edge of the cluster The lo- 
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Fig. 7. — Extinction map for cluster #91 (M82-H). The orientation follows 
that of Figs. 1 and 6 and the axes are given in arcseconds, with the x- and 
y-axis zeropoints placed roughly at the centre of the cluster The solid lines 
show isophotes as calculated from the F435W flux and the dashed circle rep- 
resents the aperture used in measuring the cluster magnitude. From this we 
deduce that the cluster is obscured by a dust lane that runs through the west- 
ern edge of the cluster. This extinction seems to be disproportionate between 
different bands, hence the shift in colour-space we see in Fig 3. 



cation and physical characteristics of this extinction patch are 
consistent with the observed shape of the cluster in the U- 
band. For cluster #91 we therefore use the age as derived by 
spectroscopy as the variable extinction should not change the 
spectroscopic features of the cluster. 

We also find tentative evidence of a similar differential ex- 
tinction pattern obscuring cluster #108, however the S/N of 
this cluster is not high enough to quantitatively assess its ef- 
fects. We are able to place the extinction range between 1.5 
and 2.0 mags in Ay in a manner that does not appear to be 
related with the underlying cluster structure. In addition, this 
cluster is not detected in the ?7-band, and this is the reason for 
its absence from the colour-colour plot (Fig 3). 

4. DISCUSSION 

In this paper we have presented an analysis of seven YMCs 
in region B of M82. Our results, along with age-dating of 
a further 29 clusters presented in Smith et al. (2007) suggest 
that these clusters are consistent with being part of a young 
cluster population that started forming <220Myr ago. We 
find extinctions in the region to vary between 1.1-2.5 mag and 
we also find differential extinction across the face of cluster 
#91 and possibly #108. These results largely contradict stud- 
ies presented in the past and we discuss them in more detail 
in the rest of this section. 

4.1. Extinction in region B 

M82 is a gas and dust rich galaxy, which is seen nearly 
edge on, causing its stellar population to be largely hidden 
in the optical. Using optic al/near-UV photometry (Smith et 
al. 2007) and spectroscopy, we find that the Ay extinction 
for clusters in our sample shows a wide range, with evidence 
of other clusters being more highly reddened. For example. 



TABLE 4 

Summary of extinction and age measurements and 
comparison to values published in 
DE Grijs et al. (2003a, denoted as G2003). 



# 


Ay A^^^""^ 

(magnitudes) 




"^spcc 


'^adopted 

(Myr) 


TG2003 


91 


2.4 


0.00 




180 


180 


5800 


97 


1.1 


0.06 


180 


190 


190 


580 


103 


1.9 


0.12 


90 


100 


100 


1000 


108 


2.5 


0.50" 




140 


140 


3200" 


113 


1.9 


1.60 


170 


100 


100 


3.8 


126 


1.5 


0.00 


260 


200 


200 


1400 


131 


1.2 


0.12 


170 


80 


80 


720 



"values taken from de Grijs et al. (2001), no data available in de Grijs et al. 
(2003a) 



cluster #108 has the highest extinction of 2.5 mag and is not 
visible in the J7-band. This variable extinction suggests that 
region B presents a line of sight into the galaxy of relatively 
low extinction, which therefore allows us to probe depths that 
are not visible on the western side of the galaxy, similar to 
a "Baade's window" in the Milky Way. Smith & Gallagher 
(2001) find a similar patch of high, variable extinction in the 
vicinity of cluster M82-F, measured by Bastian et al. (2007) 
to have Ay values between 2-4 mag, with much higher ex- 
tinction values in the surrounding region. 

The study by de Grijs et al. (2003a) concluded that the 
clusters in that sample were observable because they suffered 
from very little extinction (Ay < 1.2 mag). They suggest 
that all the clusters are situated near the surface of region B 
with respect to our line of sight, implying that the region must 
have an extremely steep gas/dust density gradient. In contrast, 
we find a larger range in extinctions and a variable overall 
extinction along this line of sight and hence on average lower 
compared to other regions of the galaxy disk. 

The selection effect created by extended areas of variable 
extinction reaching low values may explain why region B was 
thought to occupy a special place within M82, as it allows for 
the observation of a far larger number of clusters in visible 
light bands. A high number of clusters could be interpreted as 
the result of a higher cluster formation rate in the region. The 
interpretation offered in this paper is that these 'holes' in the 
dust distribution of the galaxy are in fact the source of the dis- 
crepancy in appearance between its western and eastern sides. 
This claim is supported by NICMOS near IR observations of 
the galaxy (Alonso-Herrero et al. 2001, 2003), where the two 
sides seem to have comparable luminosity and to host a sim- 
ilar number of clusters, indicating a shared cluster formation 
history. 

4.2. Cluster ages 

We find that the spectroscopic ages of the seven clusters we 
have analysed in region B are between 80 and 200 Myr Com- 
parison of the photometric and spectroscopic ages for the five 
clusters in common show that overall, the preferred spectro- 
scopic ages are towards the younger part of the allowed pho- 
tometric age range. A good example of this is cluster #126 
which has a photometric age range of 200-780 Myr (Fig. 3, 
Table 2) but the spectroscopy results clearly narrow this down 
to 160-280 Myr, and the spectrum shows that this cluster can- 
not be as old as 800 Myr 
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We now compare our ages with the ages of de Grijs et al. 
(2003a), as given in Table 4. It is clear that, with the excep- 
tion of #113, the de Grijs et al. ages are considerably older. 
In the larger photometric sample of 35 clusters presented by 
Smith et al. (2007), the ages range from 8-310 Myr** and the 
peak is at 150 Myr. Overall, we conclude that the region B 
cluster sample is much younger than presented in de Grijs et 
al. (2001,2003a). 

The ages that we find for the clusters agree well with the 
time-scale of the last encounter between M82 and M81, based 
on a value of 220 Myr (Yun 1999). de Grijs et al. (2001) 
also suggest that the age of the region B cluster population 
is consistent with this encounter but they use an earlier age 
estimate by Brouillet et al. (1991) of 5 10 Myr. 

4.3. Radial Velocities 

It has been suggested by de Grijs et al. (2003a) that region B 
is a gravitationally bound structure. In order to explain how 
such a large structure could form and survive despite the ef- 
fects of differential rotation, they suggested that the presence 
of the stellar bar (Wills et al. 2000) may have caused the for- 
mation of a circumnuclear ring, of which region B is part. It 
is well known that a ring or torus does indeed exist around 
the M82 bar, but we note that if the torus were to include re- 
gion B, it would have to have a width of '^500pc, in clear 
disagreement with the observations (r =86 pc from the centre 
according to Achtermann & Lacy 1995, WeiB et al. 2001 find 
a radius for the molecular torus of r = 65 pc). 

In order to explore the dynamical state of the region B 
cluster sample, we plot their derived radial velocities and 
the corresponding Nal interstellar absorption line mea- 
surements in Fig. 8. We have also included near in- 
frared CallA8498,8542,8662 (Ca T) stellar absorption and 
[S lll]A6531,9069 and Pa(lO) nebular emission line measure- 
ments from McKeith et al. (1993). 

Clusters #97 and 103 exhibit stellar and gas kinematics that 
are consistent with the galaxy rotation curve traced by the 
CaT velocities, and the similarity between the cluster and Na I 
measurements shows that the interstellar gas is located near 
the clusters. The radial velocities of clusters #108 and 126 
are consistent with the Ca T absorption line measurements 
and show that they are located within the main stellar disk of 
M82. However, the corresponding Nal velocities are near to 
the systemic value, indicating that a large column of cool IS 
gas is located in front of these clusters, on the outskirts of the 
disk. The radial velocities of clusters #91, 113 and 131 are 
consistent with them being located on the far- or near-side of 
the disk. The fact that the Na I IS lines are close to systemic 
argues, however, that these three clusters must be on the near- 
side of the disk otherwise we would expect the Na I velocity 
to be offset because of the greater path length. 

We thus find no indication that our sample of clusters has 
unusual kinematics; instead the large range in radial veloci- 
ties is consistent with them being located at different depths 
within the M82 disk. This also explains the large variation in 
extinction that we measure. 

Our findings allow us to test the hypothesis that region B 
is bound. In order for this to be true, the clusters should 
trace equal angles in equal amounts of time as they orbit 

* This applies to 33 of the 35 clusters in this sample. The remaining clus- 
ters #16 and 35 appear to have ages of 2.3 and 1.3 Gyr respectively, with 
minimum acceptable values of 130 and 107 Myr. These two values are poorly 
constrained and therefore not representative of the overall age distiibution. 
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Fig. 8. — Major-axis position-velocity diagram for the eastern half of 
the disk. The cluster measurements are labelled with their respective IDs. 
The CaT stellar absorption and [S III] and Pa(lO) emission measurements 
are taken from McKeith et al. (1993). 



about the galaxy centre. The galactocentric distances of 
these seven clusters places them on the gently rising part of 
the rotation curve. Assuming a constant rotational velocity 
Uiot =120 kms^^ (Fig. 8), this corresponds to an orbital pe- 
riod of 18 and 54 Myr for the inner and outer boundaries of 
the region (assumed to be marked by clusters #91 and 131 
respectively that are situated 700 pc apart). This would im- 
ply that clusters in the outermost regions should have a radial 
velocity, Wrad, three times that of the ones in the innermost 
regions. This is clearly not the case; we find no correlation 
between the cluster radial velocities and their spatial positions 
as they seem to scatter about the 120 km s^^ mark, and con- 
clude that the clusters are on regular orbits with respect to the 
stellar component of the disk. 

Finally, we remark that through a number of arguments, de 
Grijs et al. (2003a) found that the cluster population of re- 
gion B has a log-normal mass function. This supported their 
hypothesis that region B is bound, since a normal power law 
function would imply unphysically high initial densities. 

Smith et al. (2007) suggest that the turn-over in the clus- 
ter luminosity and mass functions is not real but is caused by 
the neglect of the effect of extended sources on the detection 
limit. This, in addition to the younger ages that we find, sug- 
gests that the disruption timescale derived in de Grijs et al. 
(2003a) and used in de Grijs, Parmentier, & Lamers (2005) 
is probably too short. By not imposing a short disruption 
timescale to such a young cluster population, the idea of an 
initial power-law distribution is consistent with the data. 

5. CONCLUSIONS 

We have presented Gemini GMOS spectroscopy and ACS 
t/BW photometry for seven clusters in region B of the starbust 
galaxy M82. Our aims were to obtain accurate ages, extinc- 
tions and radial velocities for the clusters. 

We have used both photometric and spectroscopic meth- 
ods to determine cluster ages. Our main age determination 
method is fitting the Balmer spectral lines with all available 
models and finding the best fit (i.e. the lowest on the over- 
all line profile fit), in order to eliminate the known Balmer 
line strength degeneracy. We find our method to agree with 
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photometric ages derived using the 3DEF method (based on 
UBVI colours). This demonstrates that the photometrically- 
derived ages are accurate but they are not as precise as the 
ones obtained from spectroscopy and confirms that the only 
truly independently reliable method for age-dating clusters in 
extragalactic environments is spectroscopy. The inclusion of 
?7-band photometry may enhance the accuracy of photomet- 
ric measurements, but unless they can be cross-checked by 
spectroscopic means, the potential age degeneracy may not 
be broken. 

We find cluster ages in the range of 80 to 200 Myr, a dis- 
tribution that is consistent with the timescale for the last en- 
counter between M82 and M81 as proposed by Yun (1999), 
namely 220 Myr. These findings in combination with the 
larger photometric sample of Smith et al. (2007), disagree 
with the 'fossil starburst' scenario proposed by de Grijs et al. 
(2001), where region B in M82 was identified with the rem- 
nant of an off-centre starburst that commenced about 600 Myr 
ago, a time scale based on their photometric determinations of 
star cluster ages. We suggest that the increased cluster forma- 
tion rate in region B is representative of the era of increased 
star formation across the galaxy disk triggered by the last en- 
counter with M8 1 . 

The extinction along our Une of sight in this extended re- 
gion appears to vary greatly, between 1.1 and 2.5 magni- 
tudes. As M82 is seen virtually edge-on, we find that M82-B 
presents a view into various depths of the body of the galaxy 
through an arrangement of 'windows' in the dust distribu- 
tion (as hypothesised for cluster M82-F by Smith & Gallagher 
2001). We also find differential extinction across the face of 
cluster #91 (and also possibly #108), which sets a serious ob- 
stacle in the photometric determination of extinction and age 
in these cases. In fact, the effect of dust in this environment is 
in some cases so grave that photometry cannot be used at all 
as an age/extinction diagnostic. 

We have also used the available spectroscopy to derive clus- 
ter kinematics. This allows us to reinforce our extinction- 
based findings and show that the seven clusters reside at dif- 
ferent depths within the disk of M82. The large scatter of clus- 
ter velocities about the gently rising component of the rotation 
curve indicates that the clusters do not move in a co-ordinated 
fashion and that region B cannot be bound. 
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Previous studies of clusters F and its neighbour L on the 
western side of the disk show that they both have an age of 
60 ± 20 Myr (Gallagher & Smith 1999; Bastian et al. 2007). 
This age fits in well with the region B cluster age distri- 
bution presented in this paper and Smith et al. (2007), and 
supports our suggestion that the increase in the cluster for- 
mation rate was not local to region B but part of a galaxy- 
wide burst. In addition, the very high masses of clusters F 
(~ IO^Mq, Smith & Gallagher 2001; Bastian et al. 2007) 
and L (4 X 10*^ Mq, McCrady & Graham 2007) imply that 
many lower mass clusters were also likely to have formed 
with F and L. 

In summary, we find region B to be optically bright owing 
to the presence of low internal extinction patches, thus 
offering a deep view into the M82 disk at radii between ~ 
400-1 200 pc. The range of cluster ages and other properties 
in this region then are typical of the evolution of the main 
body of M82 and reflect the large increase in star formation 
that occurred about 220 Myr ago when M82 last passed close 
to M81. M82 region B stands out because it is representative 
of the mid-disk zone, and is relatively clear of dust, rather 
than being a special substructure. This model also receives 
support from the larger sample of photometric ages in 
Smith et al. (2007), and will be further discussed in the 
Konstantopoulos et al. (in preparation) study of another three 
dozen M82 star clusters distributed across M82. 
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